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ABSTRACT 

We investigate the two-dimensional excitation structure of the interstellar medium in a sample of luminous infrared galaxies (LIRGs) 
and Seyferts using near-IR integral field spectroscopy. This study extends to the near infrared the well-known optical and mid-IR 
emission line diagnostics used to classify activity in galaxies. Based on the spatially resolved spectroscopy of prototypes, we identify 
in the [FeII]1.64yum/Bry - H22.12yum/Bry plane regions dominated by the different heating sources, i.e. AGNs, young main-sequence 
massive stars, and evolved stars i.e. supernovae. The ISM in LIRGs occupy a wide region in the near-IR diagnostic plane from -0.6 
to -hi.5 and from -1.2 to -hO.8 (in log units) for the [FeII]/Bry and H 2 /Bry line ratios, respectively. The corresponding median(mode) 
ratios are -hO. 18(0.16) and -h0.02(-0.04). Seyferts show on average larger values by factors ~ 2.5 and ~1.4 for the [FeII]/Bry and 
H 2 /Bry ratios, respectively. 

New areas and relations in the near-IR diagnostic plane are defined for the compact, high surface brightness regions dominated by 
AGN, young ionizing stars, and supernovae explosions, respectively. In addition to these high surface brightness regions, the diffuse 
regions affected by the AGN radiation field cover an area similar to that of Seyferts, but with high values in [FeII]/Bry that are 
not as extreme. The extended, non-AGN diffuse regions cover a wide area in the near-IR diagnostic diagram that overlaps that of 
individual excitation mechanisms (i.e. AGN, young stars, and supernovae), but with its mode value to that of the young star-forming 
clumps. This indicates that the excitation conditions of the extended, diffuse ISM are likely due to a mixture of the different ionization 
sources, weighted by their spatial distribution and relative flux contribution. The integrated line ratios in LIRGs show higher excitation 
conditions i.e. towards AGNs, than those measured by the spatially resolved spectroscopy. If this behaviour is representative, it would 
have clear consequences when classifying high-z, star-forming galaxies based on their near-infrared integrated spectra. 

Key words. Galaxies:general - Galaxies:evolution - Galaxies: structure - Galaxies:star formation - Infrared:galaxies - Infrared: ISM 
- ISM: HII regions 


1. Introduction 


The characterization of the excitation and ionization conditions 
of the interstellar medium plays a central role in our under¬ 
standing of galaxy formation and evolution across the history 
of the Universe. Galaxies are known to have different exci¬ 
tation/heating sources ranging from luminous accretion disks 
around nuclear massive black holes (AGNs), to powerful, spa¬ 
tially distributed starbursts. In general, galaxies show a complex 
excitation/ionization, two-dimensional structure due to their dif¬ 
ferent radiation field (stars versus AGNs), the presence of lo¬ 
cal (radio jets, stellar and supernovae winds) or galaxy-wide 
(tidally-induced) shocks, internal dust dumpiness, and metallic- 
ity effects (radial gradients, overall abundances). Thus, the use of 
strong emission lines for the detection and quantification of the 
ionization sources in galaxies was identified early on as an im¬ 
portant tool. The pioneering methodology of Baldwin, Phillips, 
and Terlevich (BPT; |Baldwin et al.|I98I| ), based on the use of the 


strongest optical emission lines, has been commonly used ever 
since to identify and separate AGN from star formation ([Veilleux 


& Qsterbrock|I987l|Kewley et al.|200l||Yuan et al.|20I0|) 

^axies | 


_ to de- 

fine empirical limits based on large samples of galaxies ( [Kauff- 
Imann et al. 112003 b or theoretical limits based on detailed mod¬ 


els ( Kewley et al.|2006||20I3| ). The advent of infrared satellites 
like ISO and Spitzer has also made it possible to define BPT- 
like diagnostics ([Genzel et aI.||I998||Spoon et alj2007 Groves 


et aI.|2008l|Bernard-SaIas et ai.|2009||Dale et al.|20091 Pereira- 


Santaella et aI.||2QlT] [Petrie et aL]|20IH [AIonso-Herrero et aL 


2QI2| ) based on strong mid-IR lines that are tracers of AGNs (i.e. 
high excitation lines like [NeV]I4.32yum or [OIV]25.89yum) and 
star formation (PAHs at 6.2, 7.7 and 11.3 pm). 


The study of the excitation and ionization of ISM in active 
galaxies using near-IR emission lines has been more limited. 
A BPT-like near-IR diagram involving the [Fell] 1.26yum/Pa^ 
and the H 22 . 12 /im/Br 7 line ratios was originally proposed by 
Larkin et~ar] ( |1998| ) in a study of LINERs based on low spec- 
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tral resolution (R ~ 1000) near-IR long-slit spectroscopy. Larkin 
and coworkers reported a strong linear relation in the (log-log) 
[Fell] 1.26yum/Pa^ - H22.12yum/Br7 plane, with star-forming 
galaxies showing the lowest ratios, Seyferts with intermedi¬ 
ate ratios, and LINERs with the highest ratios. Later studies 
of Seyfert 1 and Seyfert 2 galaxies ( [Rodriguez- Ardila et al. 
|2004| [2005 1 [Riffel et al] |2006| ) confirmed the trend first re¬ 
ported by Larkin and coworkers. A more recent analysis on a 
larger compilation of Seyfert 1 and 2 galaxies, including star¬ 
forming galaxies and a few LINERs ( [Riffel et ar||2013a| ), have 
shown that activity in galaxies seems to follow a unique lin¬ 
ear relation log([EeII]L26yum/PayS)= 0.749xlog(H22.12/im/Br7) 
- 0.207, with star-forming galaxies showing the lowest ratios 
and LINERs apparently the largest. 


The ionization of the ISM in luminous and ultraluminous in¬ 
frared galaxies (LIRGs and ULIRGs, respectively) shows a com¬ 
plex behaviour in the optical with the largest fraction of galax¬ 
ies classified as star forming, a large fraction (30% to 40%) as 
LINERs, and a few as AGNs ( [Veilleux et al.[[T995l [Veilleux[ 
[et al.[[199^ . The nature of LINERs in U/LIRGs has been ex¬ 
plained as the consequence of galaxy-wide merger-driven shocks 


( Monreal-Ibero et al.[[2006l [2OT0| [Rich et al.[[MTT] [MT^ , but 

other alternatives l ike a mixture of starbursts and AGN have 
also been invoked ( [Kewley et al.[[2001[ [Yuan et ar][2010[ [Kew 


[ley et al.[[2013[ ). The study of U/LIRGs in the optical includes 
additional complications derived by the amount and clumpy na¬ 
ture of their dust distribution, and by the presence of additional 
heating mechanisms due to the tidal forces associated with the 
interactions and mergers involved mostly in ULIRGs. 


A more appropriate methodology for the study of the ioniza¬ 
tion/excitation conditions of the ISM in U/LIRGs could be the 
use of near-infrared spatially resolved spectroscopy, i.e near-IR 
integral field spectroscopy (lES). This allows the investigation 
of the low- Ay, transparent ionized regions (as in the optical) 
and the more dust-enshrouded star-forming regions and/or ob¬ 
scured AGNs. These type of studies will set the ground for future 
studies with the near-IR (NIRSpec) and mid-IR (MIRI) integral 
field spectrographs onboard the James Webb Space Telescope 
(JWST). Euture observatories such as JWST will allow us to spa¬ 
tially resolve the high-redshift, infrared-luminous, star-forming 
galaxies (main sequence and above), establishing and quantify¬ 
ing the nature of their heating/ionizing sources. 


This paper presents the first spatially resolved near-infrared 
spectroscopic study of the excitation conditions of the interstel¬ 
lar medium in a sample of low-redshift LIRGs complemented 
with similar data from a sample of nearby Seyfert galaxies. Sec¬ 
tion §2 introduces the samples and near-IR integral field spec¬ 
troscopic data. The two-dimensional distribution of the ISM in 
LIRGs according to their [EellJ/Bry and H 2 /Br 7 ratios is given 
in sections §3.1, 3.2 and §3.3.1. The location of compact young 
star-forming regions, SNe- and AGN-dominated, and compos¬ 
ite starbursts-hAGN in the [EellJ/Bry - H 2 /Br 7 plane are based 
on galaxies identified as prototypes, and is presented in sec¬ 
tions §3.3.2 to §3.3.5. A discussion about possible aperture 
effects in the classification of galaxies when comparing flux- 
weighted one-dimensional with two-dimensional spatially dis¬ 
tributed classifications is given (§3.3.6). Einally, new limits and 
relations in the (log-log) [EellJ/Bry versus H 2 /Br 7 plane are 
proposed to discriminate the different excitation associated with 
young and aged stars, AGNs, and composite objects (§3.3.7). 
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Table 1: LIRG sample 


IDl 

Common 

(1) 

ID2 z 

IRAS 

(2) (3) 

Dl Scale log (Lir/Lq) 
(Mpc) (pc/") 

(4) (5) (6) 

NGC 2369 

IRAS 07160-6215 0.010807 

48.6 

230 

11.17 

NGC 3110 

IRAS 10015-0614 0.016858 

78.4 

367 

11.34 

NGC 3256 

IRAS 10257-4338 0.009354 

44.6 

212 

11.74 

ESO 320-G030 

IRAS 11506-3851 0.010781 

51.1 

242 

11.35 

IRASF 12115-4656 IRAS 12115-4657 0.018489 

84.4 

394 

11.10 

NGC 5135 

IRAS 13229-2934 0.013693 

63.5 

299 

11.33 

IRASF 17138-1017 IRAS 17138-1017 0.017335 

75.3 

353 

11.42 

IC 4687 

IRAS 18093-5744 0.017345 

75.1 

352 

11.44 

NGC 7130 

IRAS 21453-3511 0.016151 

66.3 

312 

11.34 

IC5179 

IRAS 22132-3705 0.011415 

45.6 

216 

11.12 


Notes. Col. (3): redshift from the NASA Extragalactic Database (NED). 
Cols. (4) and (5): luminosity d istance and sc ale derived using Ned 
Wright’s Cosmology Calculator ( [Wright|2006| ) given ho = 0.70, Dm = 
0.7, Dm = 0.3. Col. (6): LiR(8-1000yum) calculated from the IRAS flux 
densities / 12 , fis, feo and /iq q ([Sanders et al.[2003[ ), using the expression 
given in |Sanders & Mirabel[{[1996 . A co mplete version of this table can 
be found in |Piqueras Lopez et aL\ { 2012[ ) 


2. Samples of galaxies and observations 

2.1. Luminous infrared galaxies 

The sample of galaxies (see Table[^ consists of ten low-z LIRGs 
located at distances of ~ 45 to 85 Mpc, and covering the IR- 
luminosity range of 11.1 to 11.7 L© (logL//?). Most galaxies are 
classified as pure star forming, while only two (NGC 5135 & 
NGC 7130) are identified as hosting an AGN in their nucleus 
in addition. We obtained spectra for all galaxies with the near- 
infrared integral field spectrograph SINEONI on the VLT, using 
the seeing-limited mode, and a resolving power of about 3000 
and 4000 for the H and K bands, respectively. Eor each galaxy, 
we generated absolute calibrated maps for the H 2 2.12pm, Bry, 
and [Eell] L64yum emission lines, tracing the distribution of the 
hot molecular, warm ionized, and partially ionized gas, respec¬ 
tively. The overall field-of-view (~ 8 x 8 arcsec^) covers dis¬ 
tances of 1 to 1.5 kpc radius around the nucleus of the galaxy 
with a sampling of about 30 to 50 parsec (0.125 arcsec/spaxel), 
and a physical resolution of about 200 parsec (EWHM -0.6 arc- 
sec). Eull details of the sample observations, basic reduction, and 
calibration can be found in [Piqueras Lopez et al.[ ( [2012 ). 

We included some improvements in the reduction and cali¬ 
bration process of our SINEONI data, as explained in detail else¬ 
where (Piqueras Lopez et al. 2015, in prep.). We implemented 
different techniques to reduce the spectral and spatial noise of the 
individual observations, and to match the background of the in¬ 
dividual cubes before coadding the final data cube. After the new 
calibration of the data set, we obtained new emission line maps, 
following the same procedure as described in [Piqueras Lopez 
[et al.[ ( [20T^ . As a consequence of the new reduction process, the 
relative calibration between H- and K-band data cubes might be 
slightly different when compared with previous published mea¬ 
surements. In the worst case scenario, the differences between 
the 2012 absolute calibration and the calibration that we present 
could be as high as 20% on each band. However, on average, the 
uncertainty of the absolute flux calibration lies below < 15% on 
both bands. Therefore, we assumed a conservative 15% system¬ 
atic error for all flux measurements, added in quadrature to the 
statistical error. 

To obtain extinction-corrected flux maps, we applied the spa¬ 
tially resolved internal Ay measurements from [Piqueras Lop^ 
[et al.[ ( [20T^ . These measurements are derived from the BvyfBvd 



















































L. Colina et aL: Understanding the two-dimensional ionization structure in luminous infrared galaxies. 


line ratio and the extinction law from |Calzetti et al.| ( |2000| ), and 
are available on a spaxel-by-spaxel basis for those spaxels where 
the Br^ line is detected beyond a S/N threshold of 4 (see |Pi-| 
queras Lopez et al^|2013| for further details). In those spaxels 
where the correction is not available, the median value of the Ay 
distribution is used for each individual galaxy. Strictly speaking, 
the use of Calzetti’s law has to be applied by deriving the Ay 
for the stellar continuum first and then the extinction towards 
the ionized gas using the empirically (UV-optical) derived ratio 
Ay (continuum)/Ay (ionized gas)= 0.44. However, while extinc¬ 
tion laws differ in the ultraviolet and optical, there are similar 
extinction laws in the near infrared. For example, the differences 
(in magnitudes) between [Cardelli et al.| ( |1989| ) and Calzetti’s 
laws for the wavelengths of the Patr, Br^, and Bry emission lines 
correspond to 0.00193, 0.00658, and 0.02037xAy respectively. 

We obtained integrated measurements of the line ratios on 
different regions of interest, i.e. the peaks of the emission lines, 
[FeII]1.64yum, H 2 l-OS(l) and Bry, the main nucleus of the ob¬ 
jects, and integrated measurements over the field of view (FoV). 
The peak of each line is defined as the brightest spaxel on the 
corresponding emission line map. For consistency with our pre¬ 
vious work, we consider the nucleus of the galaxy as the bright¬ 
est spaxel of the K-band continuum map. Centred on the peaks, 
we defined a circular aperture of 200 pc radius to obtain an in¬ 
tegrated spectrum of each region. The spaxel-by-spaxel spectra 
within the aperture are de-rotated i.e. shifted to the same refer¬ 
ence frame in the spectral axis before staking (Piqueras Lopez 
et al. 2015, in prep.). This procedure increases the S/N of the 
lines by de-correlating the noise in the spectral dimension, since 
possible residuals from the sky line subtraction are no longer 
aligned. It also corrects from the artificial line broadening due 
to the intrinsic velocity field that arises from staking over wide 
apertures. After the spectra are extracted, we individually fit¬ 
ted the different emission lines (i.e. [Fell], H 2 , Bry and Br^) 
using a Gaussian profile on a spaxel-by-spaxel basis (see |Pi-| 
Iqueras Lopez et aL]|2012| for details), and obtained the corre¬ 
sponding line fluxes and Ay measurements using the Bry/Br^ 
line ratio. E xamples of spectra can be found in [Piqueras L6pez| 
jet al.jpOT^ , The [FeII]L644yum emission line flux could poten¬ 
tially be affected by the Brl2 L641/im line in absorption. How¬ 
ever, these two lines are well resolved in our R ~ 4000 spectra 
and therefore no major impact on the strong [Fell] emission is 
expected (see [Bedregal et al.||2009| for examples of the H-band 
spectra for NGC 5135). 

The final uncertainties on the line fluxes and extinction mea¬ 
surements are calculated using a bootstrap method of typically 
~ 500 simulations of each spectrum. This method consists of 
constructing a set of spectra that are statistically equivalent in 
terms of noise to the observed spectrum. We measured the noise 
of each observed spectrum as the rms of the residuals after sub¬ 
traction of the Gaussian profile at each line position. Taking this 
estimation of the noise into account, we constructed a total of 
- 500 simulated spectra where the lines are fitted again. We then 
obtained distributions of each parameter of the fitting (e.g. flux, 
line position/velocity, and cr), and calculated their corresponding 
uncertainties as the standard deviation of each distribution. The 
main advantage of this technique is that it not only takes the pho¬ 
ton noise into account, but also uncertainties due to an inaccurate 
continuum estimation or poor line fitting in low S/N spectra. 


2.2. Nearby Seyfert galaxies 

The LIRG sample has been complemented with a small sam¬ 
ple of nearby Seyfert galaxies (see Table for which high 


Table 2: Seyfert and star-forming galaxy samples 


IDl 

ID2 

z 

Dl 

Scale 

Common 

IRAS 


(Mpc) (pc/") 

(1) 

(2) 

(3) 

(4) 

(5) 

Mrk 1157 

IRAS 01306+3524 

0.0I5I67 

61.8 

291 

Mrk 1066 

IRAS 02568+3637 

0.012025 

49.2 

233 

ESO 428-G014 

IRAS 07145-2914 

0.005664 

27.1 

130 

Mrk 79 

IRAS 07388+4955 

0.022189 

98.3 

456 

NGC 4151 

IRAS Z12080+3940 

0.003319 

17.8 

86 

NGC 693 

IRAS 01479+0553 

0.005227 

18.3 

88 

UGC 2855 

IRAS 03431+6958 

0.004003 

15.9 

77 

NGC 1569W 

IRAS 04260+6444 

-0.000347 

1.45 

7 

M82 

IRAS 09517+6954 

0.000677 

4.03 

19 

NGC 4303 

IRAS 12194+0444 

0.005224 

27.5 

131 

Mrk 59 

IRAS 12566+3507 

0.002616 

14.7 

71 

NGC 7771 

IRAS 23488+1949 

0.014267 

56.6 

267 


Notes. Col. (3): redshift from the NASA Extragalactic Database (NED). 
Cols. (4) and (5): Luminosity distance and s cale derived using Ned 
Wright’s Cosmology Calculator ( |Wright|2006| ) given ho = 0.70, Dm = 
0.7, Dm — 0.3. 


spectral resolution (R - 5000) near-IR IFS data obtained with 
NIFS/Gemini are already available. The sample consists of three 


Seyfert 2 galaxies (Mrk 1066 & Mrk 1157; 

Riffel & Storchi- 

[Bergmann[ 

201 la Riffel & Storchi-Bergmann 

|2011bJ ES0428- 

G014, Riffel et al.[[2006j), one Seyfert 1.5 ga 

laxy (NGC 4151, 
2 rt 1 galaxy (Mrk 
ons of 300 to 700 

[Storchi-Bergmann et al. |2009[), and one Seyfi 
79, Riffel et al.|2013b). These data cover regi 


parsec (radius) around the AGN with an angular resolution from 
about 10 parsec (ESO428-G014 and NGC 4151) up to about 
100 parsec (Mrk 79). One of the Seyferts (Mrk 1066), while 
dominated by the AGN in the central regions, also has circumnu- 
clear star-forming knots within the central 400 pc around the nu¬ 
cleus ([Riff el & Storchi-Bergmann|2011al|Ramos Almeida et al. 
[ 2 ^ . 

The [FeII]L64yum/Bry line ratios for the Seyfert sample 
are derived from the original [FeII]/lL26)Um and Pa^ line 
emission maps. To calculate the [FeII]/lL64yum/Bry line ra¬ 
tio from the [FeII]/lL26/im/Pa(ryS, we assumed the the oretical 
ratio [FeII]/lL64yum/[FeII]/ll.26yum= 0.7646 (see e.g. [Colina 


[1993[ ), and the case B recombination factors at T=1 0,000 K 


and ne = 10"^ cm ^ ( Qsterbrock & Ferland 2006| ). There¬ 


fore, [FeII]L64yum/Bry = 4.4974 x [FeII]L26yum/Pa^. Since the 
emission lines used in the original line ratios lie very closed in 
wavelength, no correction for internal extinction effects has been 
applied to the Seyferts. 


2.3. Nearby star-forming galaxies 

There are a few nearby star-forming galaxies with published 
near-IR IFS. These are included in the analysis for a direct com¬ 
parison of the excitation properties of the brightest compact 
emitting regions in LIRGs with the nucleus and star-forming 
clumps in pure star-forming galaxies. The sample (see Table 
includes the nucleus, two off-nucleus regions (35 x 35 par¬ 
sec, each), and the central (260 x 160 parsec) region in M82 
( Forster Schreiber et al.[[200T] ), as well as one low-metallicity 
blue compact dwarf galaxy (Mrk 59, [Izotov et al.[[2009| ) and 
a small sample of nearby star-forming galaxies (NGC 1569W, 
NGC 693, UGC2855, and NGC 7771; [Dale et al.|2004D , cover¬ 
ing several morphologies and metallicities. For all these galax¬ 
ies, we obtained the [FeII]/Bry and H 2 /Bry line ratios from the 
published line fluxes for the different regions corrected by the 
specific internal extinction given in these works. 
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In addition to these galaxies from the literature, the line ratios 
for the nucleus and a circumnuclear UV-luminous star cluster 
(cluster G in |Colina et al.|2002) in the galaxy NGC 4303, located 
in the Virgo cluster, are also included. The nucleus of this galaxy, 
optically classifie d as LINER - low-luminosity Seyfert 2 ( Colina| 
& Arriba^|l999t , contains a massive (10^ Mo), UV-luminous, 


4 Myr old star cluster ( [Colina et al.||2002|), and a hard X -ray 
emitting low-luminosity AGN ( Jimenez-Bailon et al.|20Q3| ). 


3. Results 

We generated maps for the observed Bry, H22.12yum, and 
[FeII]1.64yum emission for each individual LIRG. Spatial infor¬ 
mation is typically available for a few to several thousands of 
spaxels per galaxy with a total number of over 28000 spaxels for 
all LIRGs in the sample. Since internal extinction is very clumpy 
in LIRGs and shows a wide range of values ( [Piqueras L6pez| 
et al.|2013] ), we performed subsequent analysis of the line ratios 
after applying an internal extinction correction on a per-spaxel 
basis for a large fraction of spaxels (64%) with a good determi¬ 
nation of the Ay based on the Brd/Bry ratio. For the rest of spax¬ 
els in each galaxy, the median value of the corresponding dis¬ 
tribution is applied. The results for each LIRG as well as for the 
entire sample are given in Tables[^&|^for the [Fell] L64/im/Bry 
and H 22 . 12 yum/Bry ratios, respectively. 

In addition to the maps, we also calculated the 
[FeII]L64yum/Bry and H22.12//m/Bry line ratios in each 
galaxy for a few high surface brightness regions, mainly the 
nucleus of the galaxy identified as the photometric peak in 
the SINFONI K-band line-free continuum, and the emission 
peaks of the ionized, partially-ionized, and hot molecular gas 
components as traced by the Bry, the [FeII]L64yum, and the 
H 2 2.12/im lines respectively (see Tables |7] & [^. Finally, the 
flux-weighted line ratios over the entire FOV are also given for a 
direct comparison with one-dimensional ratios for large sample 
of galaxies. 

We obtained the corresponding emission line ratios for 
the sample of Seyfert galaxies from the original NIFS emis¬ 
sion line maps (see §2.2), following the [Fell] 1.26yum/Pa^ to 
[FeII]L64yum/Bry as explained in §2.2 (see Tables [^&[^ for re¬ 
sults). 


3.1. The characterization of the ISM in LIRGs according to 
[FeiiJI.64fim/Bry ratio 

The two-dimensional distribution of the ISM in LIRGs are char¬ 
acterized by an observed (median) [FeII]L64yum/Bry ratio of 
about 0.9 that translates into an (spaxel-by-spaxel) extinction 
corrected ratio of about 1.5. Except for a galaxy (NGC 2369) 
that appears to be edge-on and heavily absorbed, the extinction- 
corrected 2D distributions appear as quasi-universal, i.e. has no 
dependence with galaxy, showing a narrow range of values with 
most regions clustering within factors less than 2 from the me¬ 
dian. There are four galaxies where only for specific regions, the 
extinction corrected [FeII]L64yum/Bry ratio is above 5 (see Ta¬ 
ble [ 7 ]). These correspond to the nuclei of NGC 7130 (Seyfert 2) 
and NGC 3110, to the luminous circumnuclear [Fell] clump in 
NGC 5135 ( [Colina et al.|2012| ), and to the overall circumnuclear 
emission in the edge-on galaxy NGC 2369. 

Starburst models, based purely on the number of ionizing 
photons and the [Fell] emission associated with the expected 
supernova rate, predict [FeII]L64yum/Bry ratios between 0.1 
and 1.4 for a continuous star formation characterized by dif¬ 
ferent IMFs ( |Colina|p^993| ). According to these models, only 
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star formation with relative steep IMFs (Salpeter or steeper) and 
low upper-mass limit (-25-30 M©), would be able to produce 
[FeII]L64/im/Bry ratios of about 1. This situation is unlikely 
and other effects not taken into account in these models like 
supernova-induced shocks and X-ray emission in aging instan¬ 
taneous bursts, or contamination by the radiation & shocks from 
an AGN, could be playing a relevant role. In fact, our Seyfert 
galaxies have typical ratios (median & mode) that are 2.6 times 
larger than those of LIRGs, indicating that the additional excita¬ 
tion mechanisms associated with AGNs, i.e. X-rays and shocks, 
could enhance the emission of [Fell] with respect to the Bry even 
in the highly ionized environments around the AGN. 


3.2. The characterization of the ISM in LIRGs according to 
the H^/Bry ratio 

The ISM in LIRGs is characterized by a (median) H 2 /Bry ratio 
of about 0.8 (see Table for details). This ratio is typically a 
factor ~ 1.8-2 lower than the corresponding ratios for Seyferts 
(see Table [^. In addition, while Seyferts cover a wide range of 
H 2 /Bry values from about 0.4 to 4.0, LIRGs share similar lower 
limits (0.3 for LIRGs) but have an upper limit of about 2.4 (95% 
of spaxels are below this value, see Table |^. This suggests that 
the excitation mechanism is such that it transitions from LIRGs 
to Seyferts, i.e. from star-forming dominated galaxies (seven out 
of nine LIRGs) to AGN-dominated galaxies. In fact, the H 2 /Bry 
values are consistent with those of Seyferts for the coronal [SiVI] 
emitting regions in the two LIRGs (NGC 5135 and NGC 1730) 
where this high-excitation AGN tracer has been detected (see 
Table 1^. For these two galaxies, the H 2 /Bry ratio in the [SiVI] 
emitting regions is about 1.5-2 times higher than in the [SiVI] 
non-emitting regions. 

Detailed studies of the near-IR H 2 line ratios in Seyferts, star- 
formi ng galaxies ([Rilfel et al.||2013a| ), LIRGs (jBedregal et al 


2009 


[Emonts et al.||2014 ), and nearby spirals ( jMazzalay et al. 


20131 indicate that the hot molecular gas is close to thermally 


excited because of a combination of X-ray radiation and shocks. 
On the other hand, the Bry line is a direct tracer of the UV- 
ionizing radiation and therefore the H 2 /Bry ratio gives an em¬ 
pirical first indication of the relative importance of each mecha¬ 
nism (UV radiation versus X-ray radiation or shocks). This will 
be discussed further in sections §3.3.2 to §3.3.4 when presenting 
the line ratios for the galaxies selected as prototypes. 


3.3. Understanding the near-IR emission line diagnostic 
diagram for star-forming galaxies and AGNs 

3.3.1. LIRGs in the [Fell]1.64yum/Bry - H 2 2.12yum/Bry plane 


The distribution of all spatially resolved regions in our sam¬ 
ple of LIRGs (~ 28000 spaxels in nine LIRGs) covers a well- 
defined region in the [FeII]/Bry - H 2 /Bry plane (see Figure [^, 
ranging from pure star-forming regions to AGNs (Seyferts), ac¬ 
cording to our empirical classification derived from the spatially 
resolved spectroscopy of bright AGN, compact young and aged 
star-forming clumps in three of our LIRGs, respectively (see sec¬ 
tions 3.3.2 to 3.3.4). 

In addition no regions simultaneously show line ratios 
log(FeII]/Bry) >1.0 and log(H 2 /Bry) > 0.8, i.e. LINERs accord¬ 
ing to the classification scheme proposed by Rilfel and collab¬ 
orators ([Rilfel et al.||2013a|). These results come as a surprise 
as all LIRGs but two (NGC 5135 and NGC 7130) in the sam¬ 
ple are classified as star forming and do not show evidence of a 


luminous AGN, according to the mid-IR (Alonso-Herrero et al. 
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Table 3: [FeII]1.64yum/Br7 line ratio for all LIRGS. Statistical results from SINFONI/IFS measurements^ 


Object 

Median 
Obs Corr 

Mode 

Obs Corr 

Mean 

Obs Corr 

(T 

Obs Corr 

P 5 

Obs Corr 

P 95 

Obs Corr 

Nspaxel 

NGC2369 

0.61 

3.58 

0.80 

4.02 

0.63 

3.99 

0.23 

2.36 

0.31 

1.40 

1.04 

7.92 

1583 

NGC3110 

0.72 

1.59 

0.69 

1.73 

0.77 

1.66 

0.23 

0.51 

0.49 

0.96 

1.18 

2.60 

2719 

NGC3256 

1.06 

1.69 

0.96 

1.82 

1.13 

2.05 

0.39 

1.46 

0.61 

0.79 

1.85 

4.64 

5553 

ES0320-G030 

0.89 

1.37 

0.91 

1.32 

0.91 

1.46 

0.22 

0.52 

0.57 

0.84 

1.30 

2.32 

4011 

IRASF12115-4656 














NGC5135 (all) 

0.95 

1.89 

1.07 

1.44 

1.21 

2.59 

0.99 

2.25 

0.49 

1.02 

3.24 

7.43 

1610 

NGC5135 ([SiVI])2 

1.06 

2.29 

1.07 

1.96 

1.40 

3.03 

1.13 

2.40 

0.48 

1.06 

4.09 

8.32 

898 

NGC5135 (no [SiVI])^ 

0.80 

1.53 

0.67 

1.37 

0.98 

2.03 

0.71 

1.92 

0.50 

0.97 

2.32 

4.93 

712 

IRASF17138-1017 

1.05 

2.00 

1.14 

2.08 

1.07 

2.08 

0.36 

0.90 

0.53 

0.83 

1.68 

3.81 

2327 

IC4687 

0.95 

1.60 

1.05 

1.76 

0.97 

1.82 

0.36 

1.01 

0.49 

0.69 

1.55 

3.87 

4212 

NGC7130 (all) 

1.18 

1.96 

0.93 

2.01 

1.42 

3.27 

0.85 

4.12 

0.50 

0.75 

3.17 

10.36 

1857 

NGC7130 ([SiVI])2 

2.08 

4.60 

2.03 

4.60 

2.20 

6.83 

0.87 

6.29 

0.89 

1.82 

3.86 

21.49 

498 

NGC7130 (no [SiVI])^ 

0.96 

1.51 

0.94 

1.24 

1.13 

1.96 

0.63 

1.53 

0.49 

0.71 

2.44 

4.77 

1359 

IC5179 

0.55 

0.82 

0.54 

0.85 

0.58 

0.94 

0.18 

0.63 

0.36 

0.42 

0.92 

1.72 

4490 

LIRGs 

0.87 

1.53 

0.95 

1.46 

0.94 

1.94 

0.49 

1.75 

0.43 

0.62 

1.72 

4.53 

28362 


Notes. ^ Columns indicate statistical results for the distribution of the observed ratios (obs), and of the extinction corrected values (corr). Last 
column (Nspaxeis gives the number of spaxels for each galaxy. ^ Statistics using only the spaxels were [SiVI] 1.96 fim emission line has been detected. 
^ Statistics using all the spaxels within the field-of-view with no [SiVI] 1.96 jjm emission line detected. 


Table 4: H 2 2.12//m/Br7 line ratio for all LIRGS. Statistical results from SINFONI/IFS measurements^ 


Object 

Median 
Obs Corr 

Mode 

Obs Corr 

Mean 

Obs Corr 

Obs 

cr 

Corr 

P: 

Obs 

Corr 

P95 

Obs Corr 

Nspaxel 

NGC2369 

1.08 

1.21 

1.18 

1.32 

1.19 

1.33 

0.59 

0.66 

0.55 

0.60 

2.30 

2.57 

1583 

NGC3110 

0.80 

0.84 

0.80 

0.80 

0.85 

0.89 

0.33 

0.34 

0.43 

0.46 

1.46 

1.53 

2719 

NGC3256 

LOO 

1.04 

1.09 

1.28 

1.20 

1.24 

0.82 

0.84 

0.33 

0.34 

2.90 

2.96 

5553 

FS0320-G030 

0.74 

0.76 

0.69 

0.69 

0.99 

1.02 

0.84 

0.86 

0.35 

0.36 

2.46 

2.53 

4011 

IRASF12115-4656 

0.69 

0.73 

0.74 

0.74 

0.73 

0.77 

0.25 

0.26 

0.39 

0.41 

1.20 

1.25 

3526 

NGC5135 (all) 

1.31 

1.37 

1.88 

1.88 

1.45 

1.52 

0.84 

0.87 

0.48 

0.50 

2.89 

3.01 

1610 

NGC5135 ([SiVI])2 

1.46 

1.53 

1.99 

1.98 

1.54 

1.61 

0.76 

0.79 

0.55 

0.58 

2.82 

2.91 

898 

NGC5135 (no [SiVI])^ 

1.09 

1.14 

0.68 

0.80 

1.34 

1.40 

0.92 

0.95 

0.44 

0.46 

2.93 

3.06 

712 

IRASF17138-1017 

0.65 

0.68 

0.58 

0.93 

0.76 

0.79 

0.45 

0.47 

0.22 

0.23 

1.61 

1.68 

2327 

IC4687 

0.540 

0.56 

0.61 

0.61 

0.57 

0.59 

0.29 

0.31 

0.17 

0.18 

1.13 

1.17 

4212 

NGC7130 (all) 

1.05 

1.10 

1.11 

1.11 

1.42 

1.47 

1.24 

1.28 

0.26 

0.27 

4.20 

4.33 

1857 

NGC7130 ([SiVI])2 

1.79 

1.88 

1.23 

1.22 

2.31 

2.42 

1.42 

1.44 

0.95 

1.06 

5.15 

5.31 

498 

NGC7130 (no [SiVI])^ 

0.76 

0.79 

0.61 

1.11 

1.09 

1.12 

0.98 

1.01 

0.23 

0.24 

3.12 

3.21 

1359 

IC5179 

0.71 

0.73 

0.84 

0.84 

0.78 

0.80 

0.37 

0.38 

0.32 

0.33 

1.46 

1.49 

4490 

LIRGs 

0.77 

0.80 

0.79 

0.67 

0.97 

1.01 

0.73 

0.75 

0.29 

0.30 

2.30 

2.38 

31888 


Notes. ^ Columns indicate statistical results for the distribution of the observed ratios (obs) and of the extinction corrected values (corr). Last 
column (Nspaxeis )gives the number of spaxels for each galaxy. ^ Statistics using only the spaxels where [SiVI] 1.96 fim emission line has been 
detected. ^ Statistics using all the spaxels within the field-of-view with no [SiVI] 1.96 jim emission line detected. 


2012| ) and X-ray emission ( [Pereira-Santaella et al.|2Qll] ). In ad¬ 
dition, according to the optical emission line diagnostics, the ex¬ 
citation of the exte nded emission is characteristic of s tar forma¬ 
tion and LINERs ( jMonreal-Ibero et J^|2006| |2010j ). This ex¬ 
tended LINER emission in LIRGs is explained by the presence 
of large scale shoc ks (jMonreal-Ibero et al.||2010| ) and compos¬ 
ite AGN-starbursts ( jKewley et al.|2001t Yun, Kewley & Sanders 
2010), and likely not due to the existence of a weak AGN. In 
the following, we further investigate these apparent inconsisten¬ 


cies by studying in more detail two LIRGs that are identified as 
prototypes of excitation dominated by star formation (IC 4687), 
by a Seyfert 2 nucleus (NGC 7130), and a third LIRG (NGC 
5135) classified as composite (AGN, and star-forming clumps 
of different ages and masses). In this galaxy, the different ex¬ 
citation mechanisms (i.e. AGN and star-forming clumps of dif¬ 
ferent ages) are at play. This will help first to spatially resolve 
and identify the dominant ionization mechanism in each region, 
and secondly define new limits and relations in the near-IR plane 
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Table 5: [FeII]1.64/im/Br7 line ratio for Seyferts. Statistical results from NIFS/IFS measurements^ 


Object 

Median 

Mode 

Mean 

cr 

P 5 

P 95 

Nspaxel 

Mrkll57 

5.43 

4.99 

7.05 

5.67 

1.43 

18.63 

868 

Mrkl066 

3.61 

3.49 

3.70 

1.16 

1.97 

5.81 

1092 

ESO428-G014 

4.15 

4.62 

4.34 

1.21 

2.83 

6.91 

60 

Mrk79 

2.80 

3.46 

4.48 

5.53 

0.73 

15.47 

1395 

NGC4151 

4.30 

4.43 

4.56 

2.15 

1.68 

8.50 

421 

Seyferts 

3.67 

3.98 

4.85 

4.57 

1.00 

13.24 

3836 


Notes. ^ Derived from the original [Fell] L26yL/m and Pa/5 ratios without further correction for internal extinction (see §2.1) 


Table 6: H 2 2.12/jmfBYy line ratio for Seyferts. Statistical results from NIFS/IFS measurements 


Object 

Median 

Mode 

Mean 

cr 

P 5 

P 95 

Nspaxel 

Mrkll57 

1.71 

1.58 

2.25 

2.07 

0.50 

5.62 

868 

Mrkl066 

1.07 

1.32 

1.13 

0.47 

0.46 

1.98 

1092 

ESO428-G014 

1.07 

1.12 

1.16 

0.40 

0.66 

2.16 

60 

Mrk79 

1.96 

2.04 

2.14 

1.10 

0.76 

4.37 

1395 

NGC4151 

0.52 

0.38 

0.78 

0.69 

0.23 

1.99 

421 

Seyferts 

1.41 

1.54 

1.71 

1.36 

0.40 

4.00 

3836 


for AGNs, pure star-forming regions (young and evolved), and 
composite objects. In our small sample of LIRGs, all galax¬ 
ies but NGC 5135 and NGC 7130 belong to the class of star¬ 
forming galaxies with no evidence of a luminous AGN. How¬ 
ever, it is known that a fraction of LIRGs have AGNs ( |Veilleux| 
|et al.|1995] ), and therefore any given LIRG will belong to one of 
the categories represented by the prototypes. 

Based on our capability of spatially resolving the distribu¬ 
tion of the emitting gas with our integral field spectroscopy, we 
differentiate five types of regions: i) bright compact Bry emitting 
regions that are identified with young (< 6 Myr), massive star¬ 
forming clumps dominated by the radiation from main-sequence 
massive stars; ii) bright compact [Fell] regions identified as aged 
(~ 8-40 Myr) star-forming clumps dominated by supernovae; iii) 
bright compact [SiVI] regions identified with nuclear AGNs; iv) 
diffuse extended [SiVI] emitting regions where the AGN radi¬ 
ation field is still detected; and v) diffuse extended non-[SiVI] 
emitting regions covering most of the galaxy. In the following 
we identify these different emitting regions in the prototypes to 
more accurately constrain the line ratios produced by different 
ionization sources. 


3.3.2. 1C 4687: prototype of LIRG dominated by star 
formation 

IC 4687 can be considered as the prototype of star-forming 
LIRG with its Bry emission (Fig.2 left upper panel) domi¬ 
nated by several circumnuclear star-forming clumps, and a weak 
Bry nucleus, which is also identified as weak Patr emitter on 
HST imaging ( [Alonso-Herrero et al.||2QQ6| . The H 2 emission 
(Fig.2 central upper panel) appears very diffuse and distributed 
over the entire circumnuclear region, with an above-average nu¬ 
cleus, and presenting the highest surface brightness. Finally, 
the [Fell] emission (Fig.2 right upper panel) follows that of 
the Bry clumpy light distribution. The overall distribution of 
the line ratios in the [FeII]L64//m/Bry - H22.12yum/Bry plane 
(Fig. lower left panel) is such that the mode of the line ra¬ 
tios is outside the previously identified region for star-forming 
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galaxies ( [Riffel et al.|2013^ but with a long tail towards low 
[FeII]L64/im/Bry (0.6) H 2 /Bry ratios (down to 0.15), character¬ 
istic of star-forming regions dominated by young, ionizing stars 
as in M82 ( [Forster Schreiber et al.|200l] see also Fig.[^. 

The bright circumnuclear Bry regions (blue points/regions 
in Fig. identify young star-forming regions with [FeII]/Bry 
and H 2 /Bry values well below 1. The nuclear region (red 
points/regions in Fig. appears as a weak and resolved ion¬ 
izing source (as tra ced by the Pacr and Bry lines), with no ev¬ 
idence of an AGN ([Alonso-Herrero et al.||2006||2012| jPereira-j 


[Santaella et al.[[MTT ), and is identified with evolved, super¬ 

nova dominated, star-forming regions (Pereira-Santaella et al. 
in prep). The line ratios are well separated from that of the 
young star-forming regions with [FeII]/Bry and H 2 /Bry ratios 
about six times higher than that of young star-forming clumps, 
respectively. The diffuse, line-emitting gas (white points/region 
in Fig.[^ has [FeII]/Bry ratios in between those of young and su¬ 
pernovae dominated star-forming regions, while the H 2 /Bry ratio 
is more consistent with that of supernovae dominated regions. 


3.3.3. NGC 7130: prototype of a LIRG dominated by a 
Seyfert 2 nucleus 

The light distribution of all emission lines in this galaxy is dom¬ 
inated by the nucleus and its surrounding regions while a few 
star-forming clumps appear along a spiral arm at distances of 
about 2.4 to 3 kpc from the nucleus (see Fig. [^. The nucleus 
itself contains a known Seyfert 2-type AGN ^Levenson et al. 


2005 [ [Alonso-Herrero et ~ar][2Q 1 2[ ) , also identified in our spec 


tra by the presence of coronal [SiVI] 1 .96pm emission line B 
[queras Lopez et aL][2012[ ). In this galaxy, there are two well- 
identified distributions in the near-IR diagnostic plane (Fig.3 
lower left panel). While a large fraction of the spatially resolved 
regions are distributed between the upper right corner of the 
star-forming region and the low AGN line ratios, there is a non- 
negligible fraction of spaxels showing very high [FeII]/Bry val¬ 
ues (> 5), as well as H 2 /Bry values similar to those of Seyferts 
(i.e. ~l-2). Note that previous optical and ultraviolet studies of 
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Fig. 1: Distribution of the spatially resolved regions (spaxel-by-spaxel) in the log([FeII]1.64yum/Br7) - log(H22.12//m/Br7) plane 
for the sample of LIRGs. Contours encircle the areas were 25%, 50%, 75%, and 90% of the spaxels are contained. The histograms 
represent the distribution for the [FeII]1.64/im/Br7 and H 2 /Br 7 ratios with mode (log) values of +0.16 and -0.17, respectively. The 
blue and green lines indicate the upper limits for the young star-forming regions and AGNs as derived from our spatially resolved 
two-dimensional spectroscopy on the high surface brightness compact regions of the prototype LIRGs (see section §3.3 for details) 


the nucleus of this galaxy have identified it as an AGN plus star¬ 
forming region. Optical emission line ratios show a smooth tran¬ 
sition from AGN-dominated in the centre to star-forming activ¬ 
ity in the external regions ( [Davies et al.||2Q14] and references), 
while the observed UV and mid-IR spectrum of the nucleus in¬ 
dicates the presence of young star clusters ( [Gonzalez Delgado] 
[et al.[[1998) and PAH emission characterist ic of star formation 
( Diaz-Santos et al.|2010[[Esquej et al.|2014[ ), respectively. 


In this galaxy, we identify four regions (Fig.3 bottom pan¬ 
els): i) the compact [SiVI] emitting region (green points/region) 
as the AGN-dominated nucleus; ii) the compact and bright Bry 
regions (blue points/regions) as the young star-forming clumps 
in the spiral arm; and iii) the diffuse AGN ionized region (yellow 
points/region) where weak [SiVI] emission around the nucleus 
is still detected; and finally iv) the diffuse ionized interstellar 
medium (white points/regions) in the extranuclear regions where 
no [SiVI] is detected, and along the spiral arm. 


The AGN-dominated nucleus clearly separates from the 
bright, compact Bry regions and from the diffuse ionized 
medium, covering exclusively the highest (i.e. > 5) [FeII]/Bry 
values. On the other hand, the bright, Bry circumnuclear regions, 
cover a range of values similar to that of the regions dominated 
by young stars in the prototype of star-forming LIRG IC 4687 
(§3.3.2). The AGN-diffuse medium clearly separates from the 
AGN- and young-star-dominated regions, covering the highest 


end of the H 2 /Bry ratio while presenting intermediate [FeII]/Bry 
values. Lastly, the diffuse non-[SiVI] emitting medium cov¬ 
ers a very wide range of H 2 /Bry values from those of com¬ 
pact young star-forming clumps to AGN-dominated nucleus, and 
with [FeII]/Bry values slightly above that for young star-forming 
regions. 


3.3.4. NGC 5135: prototype of composite LIRG 

Out of our sample of LIRGs, NGC 5135 is the only galaxy 
where clear evidence for spatially resolved and well-separated 
regions, dominated by either a Seyfert 2 nucleus, young massive 
stars, or supernova e, is well established ( [Bedregal et al. [[20091 
[Colina et al.[[2Q12[ ). Imaging of the different emission lines for 
the entire emitting region of about 1.5 kpc across, for the re¬ 
gions dominated by the AGN, young stars, supernovae, and dif¬ 
fuse ISM, and the specific integrated values for these regions are 
presented in Fig. All the regions are within a radius of 600 
parsecs from the nucleus, and are traced by the brightest cir¬ 
cumnuclear Bry (i.e. young star-forming clumps), circumnulear 
[Fell] (i.e. aged, SNe dominated star-forming clump), and com¬ 
pact [SiVI] (i.e. AGN, shown as contours in Fig. 4 upper panels) 
emitting regions. In addition, the diffuse surrounding medium 
is traced by the lower surface brightness areas outside these 
compact regions. As for NGC 7130, we can distinguish the dif- 
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Fig. 2: Maps (top panels) of the ionized (Bry, warm molecular (H 2 2.12yum) and partially ionized ([Fell] 1.64yum) gas for IC 4687, 
prototype star-forming LIRG, obtained from the SINFONI data cubes. Based on the emission line maps, the high surface brightness 
regions identified as young star-forming (in blue), aged, SNe-dominated (in red), and extended diffuse (in white) are shown for 
reference (lower right panel). The distribution of these regions in the log([FeII]1.64yum/Bry) - log(H 2 2.12yum/Bry) plane (lower 
central panel) occupy different areas indicating differences in the excitation conditions. The overall (spaxel-by-spaxel) distribution 
over the entire SINFONI FoV, as well as the line ratios for the nucleus and brightest emission line peaks are also shown (lower left 
panel). 


fuse AGN-contaminated medium traced by the presence of ex¬ 
tended, low surface brightness [SiVI], and the general diffuse 
ionized medium where no trace of coronal gas is detected. NGC 
5135 represents therefore a prototype of the so-called compos¬ 
ite galaxies where AGN, young stars, and supernovae clumps, as 
well as AGN-contaminated and general diffuse ionized medium, 
are all within about 1-2 kpc, therefore unresolved (even with AO) 
in intermediate and high-redshift galaxies. As for IC4687 and 
NGC 7130, the comparison of the distributions of the emission 
line ratios for different, spatially resolved regions (Fig. 4) give 
very important clues: 

i) The young (< 6 Myr) star-forming clumps dominated 
by massive, ionizing stars (blue points/regions) appear to have 
[FeII]/Bry and H 2 /Bry ratios factors of about six and three lower 
than that of the aged, SNe-dominated, star-forming clump (red 
points/region), respectively. These young clumps occupied a re¬ 
gion similar to that of the young stars dominated IC4687 galaxy 
and the star-forming clumps in the NGC 7130 spiral. 

ii) The aged (~ 8-40 Myr), SNe-dominated circumnuclear 
star-forming clump (red points/region) occupies a region with 
very high [FeII]/Bry ratios (~ 3 to 16). The lower range is similar 
to the ratios measured in the brightest radio supemovae remnants 
of M82 ( [Greenhouse et al.|1997| ). The upper values are close to 
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those detected in supernovae remnants, such as the Crab nebula 
( [Graham et al.||1990| ) and IC 443 ( [Graham et al.[[1987] ), where 
FeII]/Bry ratios well above 10 and up to 50 have been measured. 

in) The AGN-dominated nuclear region (green 
points/region) has an H 2 /Bry ratio similar to that of the 
SNe-dominated clump while having a [FeII]/Bry ratio ~ 3.5 
lower than that of the SNe-dominated region, i.e. close to a 
factor 2 higher than that of the young star-forming clumps. 

iv) The AGN diffuse region (yellow points/region) traced by 
the low surface brightness [SiVI] emitting gas occupies mostly 
the region of AGN-dominated ionized gas, with some overlap 
with the SNe-dominated and well disentangled from the young 
star-forming clumps. 

v) The diffuse ionized interstellar medium in the external re¬ 
gions (white points/region) has an [FeII]/Bry ratio distribution 
in between that of AGN and young star-forming clumps, and 
an H 2 /Bry distribution wider than any of the other regions but 
covering mostly the range of values characteristic of AGN- and 
SNe-dominated star-forming clumps. 

All these characteristics can be understood as variations on 
the properties of the radiation field, presence of shocks, and iron 
abundances on scales of hundred of parsecs or less. The dif¬ 
ferences between the aged (~ 8-40 Myr) SNe-dominated and 
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Fig. 3: Same as Fig. except for the AGN prototype LIRG NGC 7130. Note that the coronal [SiVI] 1.96yum emission line is 
detected in the spectra of this galaxy (contours in top panels) and therefore two new regions, AGN compact (green points) and AGN 
diffuse (yellow points) in the nucleus, and around it have been identified and occupy a region in the log([FeII] l.bdyum/Bry) - log(H 2 
2 .12/im/Br7) plane very different from that of young star-forming clumps. Also note that for NGC 7130, the nucleus coincides with 
the peak emission for all lines. 


the young (< 6 Myr) star-forming clumps can be understood as 
the combination of three different effects: i) enhancement of the 
iron abundance and emission due to supernovae (i.e. increased 
[Fell]); ii) a larger relative amount of ionizing photons in young 
clumps (i.e. increased Bry); and in) a larger relative importance 
of shocks and X-ray emission associated with supernovae (i.e. 
increased H 2 ). On the other hand, the differences between the 
AGN-dominated region and the young star-forming clumps can 
be understood as due to a larger amount of X-rays with respect to 
UV-ionizing radiation in the AGN than in young, massive star¬ 
forming clumps. An excess of X-rays penetrate more deeply in 
the surrounding medium, producing a larger partially ionized re¬ 
gion where [Fell] and H_2 emission is enhanced. Detailed stud¬ 
ies of the X-ray emission in NGC 5135 ( [Colina et al.|2012| ) es¬ 
tablished the presence of hard X-ray emitting regions associated 
with the Seyfert 2 nucleus and the bright SNe-dominated [Fell] 
emission peak, with a diffuse soft X-ray component covering the 
entire circumnuclear region. Moreover, the K-band H 2 emission 
line ratios for the nucleus and the SNe-dominated region are con¬ 
sistent with the predicted thermal ratios due to X-rays heating. 


while the ratios for young star-forming clumps depart from ther¬ 
mal ( jBedregal et al.|2Q09 1. 


The excitation properties of the diffuse gas, on the other 
hand, can be understood qualitatively as a consequence of the 
leaking of different radiation fields associated with different, spa¬ 
tially distributed ionization sources. In fact, the H 2 /Bry ratio ap¬ 
pears to be more likely associated with the leaking of X-rays 
from the AGN and SNe-dominated clump, while the [FeII]/Bry 
seems to be more associated with the UV-radiation field, as¬ 
suming the same iron abundances for the nucleus, young star¬ 
forming clumps and diffuse medium. While there is no evidence 
for radio jets in this galaxy, shocks associate d with jets can also 
enhance the [Fell] emission as in NGC 4151 ( jStorchi-Bergmann 
jet al.||2009| ), and therefore this mechanism should also be con¬ 
sidered in AGNs with energetic radio jets. 
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Fig. 4: Same as Fig.[^ except for the composite prototype LIRG NGC 5135. This galaxy is the only example where all different 
types of regions are identified: circumnuclear young, and SNE-dominated, star-forming clumps, compact and extended AGN regions 
(traced by the [SiVI] 1 .96jim emission line and shown by contours in top panels), and general diffuse ISM. All these different regions, 
i.e. excitation mechanisms, occupy different and well-separated areas in the log([FeII]1.64yum/Br7) - log(H 2 2.12yum/Br7) plane 
(lower midle panel). 


4. Discussion 


4.1. LIRGs and Seyferts. Defining new near-IR line ratio 
limits and relations to discriminate activity in galaxies 

Based on long-slit, near-IR spectra of a sample of star-forming 
galaxies and LINERs combined with previous similar da ta for 
Seyferts ( [Rodriguez-Ardila et al.|2005[ Riffel et al.|2Q13a| ) iden¬ 
tified some areas in the [EeII]1.26/im/Pa^ and H22.12/im/Br7 
plane that discriminate among the different classes of nuclear 
activity in galaxies from star formation to Seyferts and LIN¬ 
ERs. Eollowing this study, we have transformed the original 
[Eell] 1.26yum/Pa^ limits to [PeII]1.64/im/Br7 applying the cor¬ 
responding atomic and recombination parameters and correct¬ 
ing for internal extinction as explained in section §2. These re¬ 
gions, identified with broken lines in Pigs.[^&[^ are defined as 
i) star-forming galaxies characterized by log([PeII]1.64//m/Br7) 
< -rO.42 and log(H22.12/im/Br7) < -0.4; ii) Seyfert galaxies 
characterized by -hO.42 < log([FeII]1.64yum/Br7) < -hO.95 and 
-0.4 < log(H22.12yum/Br7) < -hO.78; and in) LINERs character¬ 
ized by log( [Fell] 1.64/im/Br7) > -rO.95 andlog(H22.12/im/Br7) 
> -hO.78. 

While previous work ( [Riffel et al. 2013a and references 
therein) gave the flux-weighted emission line ratios and corre¬ 
sponding limits for Seyferts, starbursts, and LINERs as above, 
we identify new limits and relations for the different excitation 
mechanism as a function of location within the galaxy (Fig.j^for 
the prototypes, and Fig.j^for all LIRGs in the sample). These in¬ 
clude those either associated with the compact AGN (green con¬ 
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tours), (circum)nuclear young and aged, SNe-dominated star¬ 
forming clumps (blue and red contours, respectively), the ex¬ 
tended AGN-diffuse medium around the nucleus (yellow con¬ 
tours), and the general diffuse interstellar medium, not directly 
associated with any of the previous regions (black contours). 
This new approach, 2D distribution versus integrated light, is 
important because many active galaxies, including high-redshift 
galaxies, do actually have an AGN and (circum)nuclear star¬ 
forming clumps of different ages. This coexistence produces a 
mixture of ionizing and excitation conditions in the ISM due to 
the relative contribution of the different radiation fields, to the 
presence of shocks induced by the AGN, radio jets, stellar winds 
and supemovae, and to the spatial distribution of the ionizing 
sources. 

We identify the new areas in the [FeII]1.64yum/Br7 - 
H22.12/im/Br7 plane (Fig.|^ occupied by the spatially separated 
and well-resolved high surface brightness regions that we asso¬ 
ciate with different dominant ionization/excitation sources in our 
prototypes (see §3.3.2 to 3.3.4): 

- SF-young (i.e. < 6 Myr, star-forming clumps): 

log([FeII]/Brr)= 0.238 -r 0.476 x log(H 2 /Brr) in the 
region limited by -0.4 < log([FeII]/Br 7 ) < -1-0.4 and -1.2 < 
log(H2/Br7) < -0.1 

- SNe-dominated (i.e. ~ 8-40 Myr, star-forming clumps): 
log([FeII]/Br 7 )= 0.705 -h 1.000 x log(H 2 /Br 7 ) in the re¬ 
gion limited by -1-0.2 < log([FeII]/Br 7 ) < -1-1.2 and -0.4 < 
log(H 2 /Br 7 ) < -rO.4 
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Fig. 5: Distribution of the line ratios in the log(FeII]1.64yum/Br7) - log(H22.12yum/Br7) plane for the various types of regions 
identified in the three prototypes (IC4687, NGC 7130, and NGC 5135). Left panel: compact high surface brightness young star¬ 
forming clumps (blue), aged, SNe-dominated clumps (red), compact nuclear AGN (green). Right panel: diffuse AGN (yellow), and 
general diffuse medium (white). The blue and green lines indicate the new upper limits for the young star-forming regions and 
Seyferts as derived from our spatially resolved two-dimensional spectroscopy on the high surface brightness compact regions. For 
comparison, the two broken lines represent the upper limits of the line ratios for the starbursts and Seyferts previously identified by 
Riffel and coworkers ( [Riffel et al.|2013a| ) based on ID spectroscopy. Contours contain 25%, 50%, 75%, and 90% of the points for 
each class. 


- AGN-compact: log([FeII]/Brr)= 1.009 - 1.312 x 

log(H 2 /Br 7 ) within the region limited by -0.3 < 

log([FeII]/Br 7 ) < -fl.5 and -0.3 < log(H 2 /Br 7 ) < -fO.9 

In addition to the classification based on the compact, high 
surface brightness regions, limits and relations for the diffuse 
medium can also be established. The AGN-diffuse region iden¬ 
tified in NGC 5135 and NGC 7130 (yellow points in Figs. [5] & 
Fig. [6|) are all within the AGN area as defined by [Riffel et S? 
P013aj ), cover an area similar to that of the AGN-compact re¬ 
gions identified above but with high values that are less extreme 
in [FellJ/Bry (no values above 10 are detected), and are char¬ 
acterized by a linear relation log([FeII]/Br 7 )= 0.353 -f 0.458 x 
log(H 2 /Br 7 ). The slope of this relation is very close to that found 
for the young star-forming clumps but clearly departs from that 
derived for a sample of star-forming galaxies, Seyferts, and LIN- 
ERs (log([FeII]/Br 7 )= 0.446 -r 0.749 x log(H 2 /Br 7 ), ob tained 
from the int egrated [FeII]1.26/im/Pa^ and H 2 /Br 7 ratios ( jRiffelj 
et al.|2013a| . A possible cause for this discrepancy could be pro¬ 
duced by the galaxies identified as LINERs in Riffels sample 
that appear showing the highest [FellJ/Bry and H 2 /Br 7 line ra¬ 
tios. LINERs are known to have very weak emission lines in the 
near-IR and therefore the subtraction of the stellar (hydrogen) 
absorption lines and continuum is essential to obtain an accurate 
measurement of the emission line fluxes as already pointed out 
previously ([Larkin et al.|1998l|Alonso-Herrero et al.|2000|[Riffel| 
[et al.|2013a[ ). The net effect would be an increase in the flux of 

the FaJS and Bry lines therefore bringing down the corresponding 
line ratios, and therefore decreasing the slope of the linear rela¬ 
tion. The regions identified as AGN compact and AGN diffuse 
agree very well with the distribution observed in Seyferts, with 
the mode values of their distributions in the near-IR diagnostic 
plane in agreement within 0.2dex (see Fig. 6 right panel). 


Finally, the diffuse interstellar medium (black contours in 
Fig. 6) has a very wide range of values covering those corre¬ 
sponding to all different excitation mechanisms. This suggests 
that the general heating and radiation field in these extended re¬ 
gions is a mixture of that of the individual heating sources (i.e. 
young stars, supernovae, AGN, shocks, etc.) weighted by their 
relative flux contribution and spatial distribution. 

This new classification scheme allows us not only to discrim¬ 
inate empirically the excitation/ionization dominated by young 
stars from that dominated by AGN but also adds a new class, 
the SNe-dominated (i.e. aged) star-forming regions. In addition 
it establishes that the excitation/ionization conditions of the ex¬ 
tended, diffuse medium is a complex mixture of all the different 
mechanisms at work, weighted by their relative contribution to 
the total energy output and also weighted by their spatial distri¬ 
bution. 


4.2. Excitation mechanisms: bright, compact line-emitting 
regions 

The line ratios for the nucleus, and the brightest Bry, [Fell], and 
H 2 regions for each LIRG in the sample have also been computed 
(Fig. left) and compared with the 2D distributions for LIRGs 
and Seyferts (blue and red contours, respectively). Note that for 
some of the galaxies, the regions emitting the brightest lines are 
spatially coincident within the available angular resolution. In 
these cases, only one of the regions is represented, typically the 
nucleus, which is also in most cases the H 2 peak region (see 
Tables [^&[^ for the specific values and details about the spatial 
coincidence of the emitting regions). In addition, line ratios for 
the nucleus and star-forming regions in nearby galaxies covering 
different metallicities and morphologies (see §2.3 for details) are 
also included in this analysis for a direct comparison. 


Article number, page 11 of 18 
























































A&A proofs: manuscript no. aa-SINFONI-BPT_arxiv 


DO 





- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 

log(H2/BrY) 




- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 

loglH^/Bry) 


- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 

log(H2/BrY) 


Fig. 6: Distribution of the line ratios for all LIRGs in the log(FeII]1.64/im/Br7) - log(H22.12yum/Br7) plane for different types of 
regions: high surface brightness clumps and nucleus (left panel), for the diffuse, AGN and non-AGN, extended regions (central 
panel), and for the compact and diffuse AGN medium in LIRGs and low-z Seyferts (right panel). The blue and green lines indicate 
new upper limits for young star-forming regions and Seyferts as derived from the compact, high surface brightness regions in our 
LIRG prototypes (see Fig.[^. For comparison, the two broken lines represent the limits for the starbursts and Seyferts previously 
identified by Riffel and coworkers ( [Riffel et al.|2Q13a| ) based on ID spectroscopy. Contours contain 25%, 50%, 75%, and 90% of 
the points for each class. 


Most of the brightest emitting regions are located either close 
to the mode of the 2D distribution of the emission line ratios 
for LIRGs (blue contours), or towards lower [FellJ/Bry and 
H 2 /Br 7 ratios. These agree with the range of values observed 
in the nucleus and off-nuclear regions of M82, in the circum- 
nuclear young star cluster in NGC 4303, and in the nucleus of 
star-forming galaxies with metallicities around solar, and are 
located well above the values derived in low-metallicity (Zq/3 
or less), star-forming galaxies (Mrk 59 and NGC 1569W). A 
small fraction of regions are located in the area previously iden¬ 
tified as SNe-dominated (nucleus of NGC 3110, [Fell] peaks in 
NGC5135 and IC 4687) and AGN-dominated (NGC 7130). The 
two galaxies with Seyfert 2 nucleus show a very different be¬ 
haviour. While the Seyfert 2 nucleus of NGC 7130 lies close 
to the locus of the nearby Seyferts (red contours) and close to 
the nucleus of the NGC 4303 (low-luminosity AGN plus nu¬ 
clear young cluster, see references in §2.3), the nucleus of NGC 
5135 shows a significantly lower [FellJ/Bry ratio. In addition to 
an AGN, NGC 7130, and NGC 5135 have strong starbursts in 
their central 100 p c region ([Gonzalez Delgado et al.|199~8l|Diaz-| 
Santos et al.|2010[ |Esquej et al.|2014| ). Therefore changes in the 

emission line ratios could be expected if the nuclear radiation 
field is a mixture of AGN and young stars with different relative 
contributions, and/or due to aging effects in the nuclear star clus¬ 
ters. Finally, only the nuclear regions in ES0320-G030 appear 
as outliers with very high values (~ 6) for the H 2 /Br 7 ratio, and 
intermediate [FellJ/Bry (Fig. see also Fig. [nj. We identified 
this region as H 2 -dominated requiring detailed investigation. 


In summary, most compact bright regions in LIRGs ap¬ 
pear to be similar to solar-metallicity, star-forming clumps in 
nearby galaxies, and well separated from low-metallicity, star¬ 
forming clumps. A smaller fraction of regions appear to be con¬ 
sistent with Seyfert nuclei, although the contribution of young 
stars and SNe appears to play a relevant role. The nucleus of 
ES0320-G030 is identified as an outlier to the general be¬ 
haviour because of its strong H 2 emission relative to that of Bry. 
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4.3. Spatially resolved versus integrated line ratios. 

Implications for the classification of high-z star-forming 
galaxies 

While all the analysis we performed so far is based on the two- 
dimensional distribution of the line ratios, most studies of the 
excitation mechanisms in large galaxy surveys have been based 
on the nuclear/integrated line ratios. The main differences be¬ 
tween these two measurements is that, while the integrated val¬ 
ues represent a flux-weighted measurement and therefore it is 
dominated by the brightest observed regions, the median of the 
distribution of the spatially resolved regions gives a value more 
representative of the general excitation and ionization conditions 
of the extended diffuse medium in galaxies. This is particularly 
important for LIRGs, where different ionizing sources coexists, 
and where the clumpy structure of the dust extinction can play a 
relevant role. 

A comparison of the flux-weighted (i.e. integrated) and the 
2D distribution (median) values for the LIRGs and Seyferts is 
presented (Figs. [ 7 ] right panel). For the Seyferts in the sam¬ 
ple, the integrated and median values for the line ratios agree 
within less than 0.1 dex, and therefore occupy the same loca¬ 
tion in the [FeII]/Bry - H 2 /Bry plane, near the locus of the 
2D distribution. On the other hand, clear differences appear 
for LIRGs. While IC5179 shows integrated and median values 
within O.ldex, IC4687 and IRASF17138-1017 show integrated 
values slightly lower (factors 0.5 to 0.6) than their median ratios. 
The rest of the LIRGs show larger differences in the [FeII]/Bry 
ratio with integrated over median values ~ 1.5 - 2.2 times higher. 
The net effect for these galaxies is that they appear with higher 
excitation, i.e. moving towards the locus of the AGNs, when 
the flux-weighted integrated values are considered instead of the 
median. Whether this is a general behaviour in LIRGs or galax¬ 
ies in general (similar changes in the line ratios when compar¬ 
ing nuclear and off-nucle ar regions in star-form ing galaxies have 
been reported recently by [Martins et al.|20T3] ), needs further ex¬ 
ploration with a larger volume of galaxies with available inte¬ 
gral field spectroscopy. This would have obvious implications 
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Fig. 7: Near-IR emission line ratios for the nucleus and brightest Bry, H 2 , and [Fell] emitting regions in LIRGs (left panel). Note 
that in some LIRGs these regions coincide (see Tables |^&[^ for details). The distribution of spatially resolved regions (spaxel-by- 
spaxel) are represented by contours for both LIRGs (blue) and Seyferts (red). Contours for LIRGs and Seyferts contain 25%, 50%, 
75%, and 90% of the points. Examples of nearby star-forming galaxies like M82, Mrk 59, composite AGN-fstars (NGC 4303), and 
other low-z SFG are also indicated. The emission line ratios for the sample of LIRGs and Seyferts as derived from the integrated 
values (i.e. ID flux-weighted), and median of the 2D distribution (i.e. spaxel-by-spaxel) of values is also presented (right panel). 


when classifying high-z galaxies for which the linear resolution 
would be in general very limited (~ 1.5-2 kpc), even with sub- 
arcsec angular resolution with current and future facilities like 
the JWST. 


5. Conclusions 

We investigated the two-dimensional excitation structure of the 
interstellar medium in a sample of low-z LIRGs and Seyferts, 
using near-IR integral field spectroscopy. Based on the spatially 
resolved spectroscopy of three (selected) LIRG prototypes, we 
identified in the [Fell] L64yum/Bry - H 22 . 12 yum/Bry plane the 
regions dominated by the different heating sources, including 
young main-sequence massive stars, supernovae, and AGNs. 
This extends to the near infrared the well-known optical emis¬ 
sion line diagnostics used to classify the activity in galaxies. The 
main conclusions are: 

1. The two-dimensional (i.e. spaxel-by-spaxel) distribution of 
the excitation conditions of the ISM in LIRGs occupy a 
wide region in the near-IR diagnostic plane from -0.6 to 
-1-1.5 and from -1.2 to - 1 -O .8 (in log units) for the [FeII]/Bry 
and H 2 /Bry line ratios, respectively. The corresponding me- 
dian(mode) ratios are -f0.18(0.16) and -f0.02(-0.04), respec¬ 
tively. 

2. The two-dimensional distributions of the line ratios of LIRGs 
and nearby Seyferts partially overlap, with Seyferts showing 
on average larger values by factors ~ 2.5 and -1.7 for the 
[FeII]/Br 7 (4.9 vs. 1.9) and H 2 /Br 7 (1.7 vs 1.0) ratios, re¬ 
spectively. 

3. We identified a range of line ratios for the different ioniza¬ 
tion/heating sources thanks to our integral field spectroscopy 
of well-resolved and spatially separated, high surface bright¬ 
ness, line emitting clumps: i) AGN dominated, ii) young 
stars dominated, and in) supemovae dominated. These three 
heating mechanisms occupy mostly well-separated regions 
in the log([FeII]/Br 7 ) - log(H 2 /Br 7 ) plane. 


4. New areas in the near-IR diagnostic emission line plane 
are defined for different heating sources with the follow¬ 
ing limits in the log([FeII]/Br 7 ) and log(H 2 /Br 7 ) ratios, re¬ 
spectively: AGN-dominated: [-0.3, -fL5] and [-0.3, +0.9]; 
young stars-dominated: [-0.4, +0.4] and [-1.2, -0.1]; and 
SNe-dominated: [+0.2, +1.2] and [-0.4, +0.4]. Within these 
limits, the 2D distribution of the regions dominated by 
each of the mechanisms appears to cluster around lin¬ 
ear relations given by: i) AGN-dominated: log([FeII]/Br 7 )= 
1.009 - 1.312 X log(H 2 /Br 7 ); ii) young stars-dominated: 
log([FeII]/Brr)= 0.238 + 0.476 x log(H 2 /Brr); and iii) SNe- 
dominated: log([FeII]/Br 7 )= 0.705 + 1.000 x log(H 2 /Br 7 ). 

5. We also identified the excitation conditions of the diffuse, 
extended ISM outside the bright clumps. The diffuse ISM 
occupies a large region that completely overlaps with that of 
the individual excitation mechanisms (i.e. AGN, young stars, 
and supernovae), indicating that its radiation/heating field is 
due to a mixture of the different ionization sources, weighted 
by their relative flux contribution and spatial distribution. 

6 . Compared with the median of the two-dimensional distribu¬ 
tion, the integrated line ratios in LIRGs appear compatible 
with higher excitation conditions, i.e. towards AGNs. If this 
behaviour is general, it would have clear consequences when 
classifying high-z, star-forming galaxies based on their near- 
infrared integrated spectra. A similar situation could also 
happen with the optical emission line ratios that would de¬ 
serve detailed exploration with current and future observato¬ 
ries like the James Webb Space Telescope. 
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Table 7: [Fell] 1.64yum/Br7 line ratios for LIRGs. Values for emission line peaks and integrated measurements 


A&A proofs: manuscript no. aa-SINFONI-BPT_arxiv 


o 

U 


(D 

Ph 


rO 

O 


o 

U 


Ph 


o 


o 

U 


<D 

Ph 


O 


o 

U 


O 


o 

U 


X 

O 


in 

p 

oo 

o^ 

X 

oo 

o 

o 

cn 

cn 

G) 

p 


o 

G\ 

(N 

ci 


d 

cn 

d 

d 

cn 

d 

+1 

+1 

i +1 

• +1 

+1 

+1 

+1 

+1 

(N 

X 



oo 

X 

p- 

oo 



p 

p 

G\ 

X 

p 

p 

p- 

p- 

cn 


d 

d 

in 


o 

o^ 

p- 

in 


(N 


in 


o 

in 


o 

O 


o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

: +' 

+1 

+1 

+1 

+1 

X 

in 

o 

p- 



p- 

oo 



p 

p 

X 

p- 

p- 

X 


d 


cn 

d 

d 

(N 

d 

cn 


o 


in 

(N 


(N 

G\ 

(N 

p 

X 

oo 


p 

(N 

d 

X 

(N 

d 

d 

d 

p- 

d 

+1 

+1 

i +1 

• +1 

+1 

+1 

+1 

+1 

o^ 

p- 

X 

X 

(N 

oo 

O 

p- 

p-^ 


G) 

oo 

G) 

p 

p 


(N 

X 




(N 

in 


X 

o^ 

OO 

X 

cn 

OO 


in 

o 

o 

X 

o 

o 

o 


o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

■ +1 

+1 

+1 

+1 

+1 

p- 

X 


cn 

cn 

in 

oo 

o 

oo 

fr^ 

p 

p 

p 

p 

p- 


d 

d 





(N 

d 

o^ 


cn 

p- 





(X3 

cn 

in 

cn 


o 

oq 

(N 

d 

d 

d 

d 

d 

d 

cn 

d 

+1 

+1 

: +1 

: +1 

+1 

+1 

+1 

+1 

<N 

o^ 


cn 

o^ 



o 

p- 

o^ 

p 

p 

X 

X 



(N 

d 



d 

d 

in 


X 

in 

X 

cn 


<N 


in 

O 

o 

o 

o 

o 

O 


o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +' 

: +' 

+1 

+1 

+1 

+1 

cn 

X 

p- 

p- 

in 

cn 

o 

oo 

(X3 

in 

r-- 


in 

p- 

p 

X 

d 

d 

d 

d 

d 

d 

(N 

d 

oo 


(N 


in 

cn 

O 

cn 

o^ 

a\ 

p 

X 

oo 

X 

G\ 

(N 

d 

in 

(N 

d 

d 

d 

cn 

d 

+1 

+1 

: +1 

■ +1 

+1 

+1 

+1 

+1 

o^ 

o 

(N 

X 

(N 


p- 

o 

p-^ 

G) 

p 

oo 

G) 

p 


p 

(N 

in 

cn 



(N 

in 


X) 

o^ 

o 

X 

cn 

X 


in 

o 

o 

X 

o 

o 

o 


o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +' 

• +1 

+1 

+1 

+1 

+1 


X 

p- 

cn 

cn 

cn 

p- 

o 

oo 

fr^ 

p 

p 

p 

G\ 

p- 


d 

d 




d 

(N 

d 


oo 



(N 

oo 

in 

oo 

<N 

p 

p 

p 

p- 


X 

cn 

in 

d 

d 

d 

d 

d 

cn 

d 

+1 

+1 

: +1 

: +1 

+1 

+1 

+1 

+1 

cn 

p- 

X 

o^ 

oo 

cn 



cn 

p-_ 

p 

p 

p 

p 

p 

o^ 

in 

cn 

ri 

(N 



cn 

d 

(N 

o^ 


O 

p- 

p- 

o^ 

X 

(N 

o 



o 

o 

o 

o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

■ +1 

+1 

+1 

+1 

+1 

X 

X 

in 

o 

(N 

in 


(N 

<N 


o^ 


o^ 

oo 

p 

X 


d 

d 


d 

d 


d 


o^ o 

\o ^ 
cn ^ 
(N cn 
U U 

o o 


X o 

CM V 
cn o 
U ^ 

O s 

^ o 

^ CO 

W 


o 

(N __I <N CO 

^ ^ ^ of 

cnoot^O.-r 
^ cn oo cn On 

IT) ^ VO ^ 


a 

(L> 

Ph 

(N 

X 

X 

p 

p 


p 


p 

'o 

p 

'o 

o 

u 

<D 

Ph 


"d 

p 

p 


p 


p 

d 

> 

'So 

05 

<D 

P 

> 

O 


S 

o 

o 

p 

_o 

o 

p .p 
•p p 

X P 
<D P 
-d 'o 

S -S 

05 


p 

o 






o _ 

p ^ 

(Z3 ' ' 

OJ P 

^ d 

P l-H 


p p 
O P 


p p 

'o 
p o 
P (D 
<D H 

d M 
p 

p d 

p _ 

d CT 


£r 


Ph ^ 
Ci d 
. "d 

!/! -p 

^ P 

^ o 

Z d 


P 

(D 

a 

d 


p 

(D 

a 

"O 

B 

p 

(PO 

d 


"O 

P 

p 

c/5 

d 

Ph 

d 

_P 

P 

o 


d 

P 

I 


o 

HH 

hJ 


p 


d 

P 


ca 




00 

d 

i 


d 


o 

U 


X 

O 


o 

U 


p 

(N 

X 


d 

P 

X 


cn 

On 

(N 

p- 

cn 

cn 

p 

p 

p- 

o 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

: +1 

+1 

p- 

On 

• cn 

• o 


p 

r-- 

p 

p 

(N 


d 

in 


d 

p- 

p- 

X 

oo 


o 

o 

cn 

o 

o 

d 

d 

d 

d 

d 

+1 

+1 

■ +1 

i +1 

+1 

in 

o 

p- 

(N 

X 

G\ 


p 

X 

(N 

d 

d 

p- 


d 


X 

O 


o 

U 


X 

O 


o 

U 


X 

O 


o 

U 


X) 

O 


o o o 
+1 +1 +1 
»n (N oo 

oH ^ ^ 


o^ 

cn 

(N 

o 


in 

X 

(N 


in 

G) 

p 

p- 




d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

+1 

+1 

+1 

+1 

+1 

oo 

(N 

X 

(N 

X 

(N 


p- 


OO 

(N 

G\ 

p 

in 


p-^ 

d 

d 

X 



d 

d 


p- 

p- 

in 

o^ 

X 

(N 

cn 

p- 

o 

o 

in 


o 

O 

o 

o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

• +1 

+1 

+1 

+1 

+1 

+1 

cn 


o 

o^ 

o 

On 

oo 




G) 

X 

p 

p- 

X 

p 

d 

d 

in 



d 

d 


oo 

cn 



o 

(N 


oo 

T-1 

T-1 

T-1 

X 

T-1 

O 

o 

X 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

i +1 

+1 

+1 

+1 

+1 

+1 


in 

On 

oo 

(N 


p- 

o 


X 

p- 

oq 

X 



p-^ 

d 

d 

d 


d 

d 

d 


p- 

(N 

cn 

o^ 

(N 



p- 

o 

O 

o 


O 

o 

o 

o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

+1 

+1 

+1 

+1 

+1 

(N 

(N 

X 

X 

o 

X 

cn 

p- 


X 

p- 

X 

X 



p 

d 

d 

d 


d 

d 

d 


o 

cn 


o 


in 

(N 

o^ 

(N 

in 

p 

p 

p- 



X 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

+1 

+1 

+1 

+1 

+1 


(N 


(N 

X 

(N 

o^ 

(N 


OO 

oq 

G) 

p 

in 

in 

P-^ 

d 

d 

in 



d 

d 


p- 

p- 

X 

o^ 

X 

(N 

(N 

p- 

o 

o 

p- 


o 

O 

O 

O 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

■ +1 

+1 

+1 

+1 

+1 

+1 

(N 

(N 


o^ 

o 

On 

X 

in 



p-^ 

X 

p 

p- 

in 

p 

d 

d 

in 



d 

d 


oo 

cn 

oo 

p- 

o 

oo 

p- 

p- 


(N 

(N 

oo 

p- 

o 

o 

o^ 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

+1 

+1 

+1 

+1 

+1 


o 

cn 

in 


o^ 

o^ 


p 

G\ 

G\ 

o^ 

p 

p- 

cn 

p 


d 

d 

d 


d 

d 


p- 

p- 

in 

oo 


(N 




o 

o 

o 

o 

O 

o 

o 

d 

d 

d 

d 

d 

d 

d 

d 

+1 

+1 

: +1 

+1 

+1 

+1 

+1 

+1 

o 

(N 

oo 

p- 

o^ 


oo 


P 

OO 

oo 

oo 

p 

p- 

cn 

p-^ 


d 

d 

d 


d 

d 



P ^ 

o^ o in 

X) ^ (N 


Fi n c ° 

U U rh cn 
0 0^0 
X X ^ 


o 

N ^ CN CO 

in I Tt ■ —L of 

cn oo o p 

^ cn oo cn On 

in ^ O) ^ 

Hpuo^ 

< z 

X 


'd 


p 

d 

"d 

'd 

p 

'o 


X 

X 


p 

d 

> 

’Ei) 

p 

c/3 

d 

13 

> 

U 

o 


"d 


g § 

^.S 

§ ° 

^ d 

X ^ 
O v/i 

> Sh 

p rTP 

d td 

_n 

O XL 

X 

c/3 P 
CD p 

d H 
p 

O P 

p s 

Ejq g 

d .p 
p o 

p d 

° a 
g 

S OJ 


P d 
d u. 

"d ^ 
d § 

p c/3 

tJD P 


d 


d 


.a a 

s § 

^:s 

-I 

oi o 

X d 


c/3 


P 


1^ 
^ Ph 


Article number, page 14 of 18 













L. Colina et aL: Understanding the two-dimensional ionization structure in luminous infrared galaxies. 


6. Appendix 

The appendix presents in graphical form (Figs. to the 
emission line maps and distribution of the emission line ratios 
in the log(FeII]1.64yum/Br7) - log(H22.12yum/Br7) plane for the 
sample of LIRGs, excluding the prototypes that have already 
been discussed (Fig s. [^to[4|). Notes on individu al galaxies can 
be found elsewhere ( [Piqueras Lopez et al.|2012| ) 
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Fig. 8: Emission line maps and near-IR diagnostic planes for NGC2369. Colour-coded symbols as in Fig.j^ 
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Fig. 9: Emission line maps and near-IR diagnostic planes for NGC3110. Colour-coded symbols as in Fig.j^ 
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Fig. 10: Emission line maps and near-IR diagnostic planes for NGC3256. Colour-coded symbols as in Fig.j^ 
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Fig. 11: Emission line maps and near-IR diagnostic planes for ES0320-G030. Colour-coded symbols as in Fig.j^ 
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Fig. 12: Emission line maps and near-IR diagnostic planes for IRASF17138-1017. Colour-coded symbols as in Fig.| 
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Fig. 13: Emission line maps and near-IR diagnostic planes for ICS 179. Colour-coded symbols as in Fig.j^ 
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